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Abstract An increasing number of experiments at the Belle, BNL, CERN, DAONE 
and SLAC accelerators are confirming the violation of time reversal invariance (T). The 
violation signifies a fundamental asymmetry between the past and future and calls for a 
major shift in the way we think about time. Here we show that processes which violate 
T symmetry induce destructive interference between different paths that the universe 
can take through time. The interference eliminates all paths except for two that repre- 
sent continuously forwards and continuously backwards time evolution. Evidence from 
the accelerator experiments indicates which path the universe is effectively following. 
This work may provide fresh insight into the long-standing problem of modeling the 
dynamics of T violation processes. It suggests that T violation has previously unknown, 
large-scale physical effects and that these effects underlie the origin of the unidirection- 
iy-j \ ality of time. It may have implications for the Wheeler-DeWitt equation of canonical 
"vi" ; quantum gravity. Finally it provides a view of the quantum nature of time itself. 
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' 1 Introduction 



The physical nature of time has been an enigma for centuries. The main tools for 
discussing its unidirectionality are the phenomenological arrows of time. A great deal 
of progress has been made in recent decades in linking various arrows together [TJ. A 
notable exception is the matter-antimatter arrow which arises from the violation of 
charge-parity conjugation invariance (CP) in meson decay CP violation, which 

was first discovered by Cronin, Fitch and coworkers [5] in 1964 in the decay of neutral 
kaons (K mesons), provides a clue to the origin of the large-scale matter-antimatter 
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imbalance of the universe [6]. The combined CPT invariance has been confirmed in kaon 
decay to a high degree of precision [7] suggesting that neutral kaon decay also violates 
T invariance. To date no violation of CPT has been observed in B meson decay [8] and 
so the CP violation in B meson decay also appears to be consistent with T violation. 
But more importantly, experimental evidence has confirmed the direct T violation in K 
meson decay [H] independently of CPT invariance. Nonetheless, T violating processes 
are relatively rare and the magnitudes of the violations are relatively small, and so 
despite their important asymmetric temporal nature, they are often regarded as having 
little impact on the nature of time. As a result, the physical significance of T violation 
has remained obscure. 

The CP violation due to the weak interaction is described in the Standard Model 
of particle physics by the Cabibbo-Kobayashi-Maskawa (CKM) matrix [TUltTT]. Never- 
theless, a fundamental and intriguing question remains open. In principle, the CKM 
model provides the basis for the time evolution of mesons in terms of Schrodinger's 
equation and an associated Hamiltonian H . The evidence of T violation implies that 
THT^ 1 ^ H where T represents the time reversal operation [12]. But how should one 
incorporate the two Hamiltonians, H for time evolution towards our future and THT^ 1 
for evolution towards our past, in one quantum equation of motion? At issue is the fact 
that Schrodinger's equation describes consistent unitary evolution in both directions 
of time: the unitary evolution by negative time intervals merely "back tracks" exactly 
the same time evolution described by positive time intervals. Yet, paradoxically, the 
evidence of T violation processes appears to suggests that Nature does not work in this 
way in general. 

Currently, to model T violating processes we must choose the direction of time 
evolution and a specific version, say H, of the Hamiltonian in order to write down 
the Schrodinger equation. By applying the time reversal operation we then obtain the 
Schrodinger equation involving the THT~ l version of the Hamiltonian for evolution in 
the opposite direction of time. But at any point we are restricted to having a dynamical 
equation of motion for a specific direction of time evolution and a corresponding specific 
version of the Hamiltonian. We do not have a dynamical equation of motion for the 
situation where the direction of time evolution cannot be specified. In this case there is 
no argument for favoring one version of the Hamiltonian over the other and so if one 
version of the Hamiltonian is to appear in the equation of motion then so must the 
other. This problem becomes critical when we attempt to describe the universe as a 
closed system as needed in cosmology, for then the direction of time evolution cannot 
be specified because a closed system precludes any external clock-like device for use 
as a reference for the direction of time. There is, therefore, no satisfactory quantum 
formalism for describing a universe that exhibits T violation processes. This failure of 
standard quantum theory is all the more serious given that CP and, by implication, T 
violation play a critical role in baryogenesis in the early history of the universe [6lfT3] . 
It is clearly an important problem that lies at the very foundations of quantum theory. 
Its resolution calls for a major shift in the way we think about time and dynamical 
equations of motion. 

We address this issue by deriving a single quantum dynamical equation which incor- 
porates the two versions of a T-violating Hamiltonian. We begin by using Feynman's 
sum over histories [TJ] method to construct the set of all possible paths that the uni- 
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verse can take through time. The set represents every possible evolution zigzagging 
forwards and backwards through time. We then show how T violation induces destruc- 
tive interference that eliminates most of the possible paths that the universe can follow. 
Effectively only two main paths, representing continuously forwards and continuously 
backwards time evolution, survive the interference. Physical evidence from accelerator 
experiments allows us to to distinguish between the two directions of time and determine 
which path the universe is effectively following. In short, T violation is shown to have 
a previously unknown large-scale physical effect that underlies the unidirectionality of 
time. This is the main result of this work. We end with a discussion that includes the 
potential impact for cosmology and the Wheeler-DeWitt equation in canonical quantum 
gravity. 

2 Possible paths through time 

Consider a system, which we shall refer to simply as the universe, whose composition 
in terms of size, matter and fields is consistent with the visible portion of the physical 
universe. For clarity we assume that the universe is closed in the sense that it does not 
interact with any other physical system and further, that there is no physical system 
external to the universe. This means that there is no clock external to the universe and 
so our analysis needs to be unbiased with respect to the direction of time evolution. 
Nevertheless it is convenient to differentiate two directions of time as "forward" and 
"backward" associated with evolution in the positive— t and negative— t directions of 
the time axis, respectively, although we stress that neither direction necessarily has 
any connection to our actual everyday experience. Also, we begin our analysis with the 
universe in a state \ipo) which we shall call the "origin state" without reference to the 
direction of time. With this in mind we express the evolution of the universe in the 
forwards direction over the time interval r as \i/jf(t)) = Up(t)\iPq) where 



is the forwards time evolution operator, and Hp is the Hamiltonian for forwards time 
evolution. Throughout this article we use units in which H — 1. If time reversal sym- 
metry were obeyed, the evolution in the opposite direction of time would be given by 
exp(iTHp)- To accommodate the time asymmetry of T violation processes, we need to 
replace Hp in this expression with its time reversed form. Thus we define the state of 
the universe after backwards evolution over a time interval of the same magnitude by 
\^b{t)) = U b {t)\iI)q) where 



and Hb = THpT~ l is the Hamiltonian for backwards time evolution. This gives 
Ub( t ) — TUp^T" 1 , in accord with Wigner's definition of time reversal |T2]. In 
other words, Hp and Hb are the generators of time translations in the positive— t 
and negative— t directions, respectively. 

The parameter r in these expressions represents a time interval that could be ob- 
served using clock devices which are internal to the universe. A simple model of a clock 
for our purposes is a device which evolves unitarily through a sequence of orthogonal 
states at sufficient regularity to give a measure of a duration of time to any desired 




(1) 




(2) 
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accuracy. The orthogonal states would represent the eigenstates of a pointer observable 
which acts a reference for time. A harmonic oscillator is a device of this kind [15]. In 
order to avoid ambiguities in the definition of time intervals for different time directions, 
we assume that the Hamiltonian _£/"( clk ) describing our clock model is T invariant and 
so T i7( clk ) T^ 1 = jy( clk ) during its normal operation. This gives an unambiguous oper- 
ational meaning of the parameter r as a time interval. We shall refer to the readings of 
such clocks as "clock time" . 

The matrix elements (4>\U f (t)\^o) and (4>\U b {t)\tPq) represent the probability am- 
plitudes for the universe in state |^o) to evolve over the time interval r to \<p) via two 
paths in time corresponding to the forwards and backwards directions, respectively. 
Given that we have no basis for favoring one path over the other, we follow Feynman 
[H] and attribute an equal statistical weighting to each. Thus, the total probability am- 
plitude for the universe to evolve from one given state to another is proportional to 
the sum of the probability amplitudes for all possible paths through time between the 
two states. In the current situation we have only two possible paths and so the total 
amplitude is proportional to 

(<P\Uf(tM ) + (0|£/ b (t)|Vo> = Wf{t) + U B (rMo) ■ (3) 

This result holds for all states \<f>) of the universe and so the unbiased time evolution of 
|^o) over the time interval r can be written as 

\#{t)) = [U f {t) + U b {t)]\1> ) (4) 

which we call the symmetric time evolution of the universe. To reduce unnecessary 
detail we use un-normalized vectors to represent states of the universe. We could have 
included in an additional phase factor e 10 to represent an arbitrary relative phase 
between the two paths, such as ((()\U f {t)\tPq) + e l6 \<t>\V "b(t)|^o) > but the factor has no 
net physical effect for our analysis and so we omit it. 

It follows that the symmetric time evolution of the universe in state |^(t)) over an 
additional time interval of r is given by 

\V(2t)) = [U f (t) + U b (t)]\V(t)) 
= [U F (r) + U B {r)f |Vo> • 

Repeating this for iV such time intervals yields 

\nNr)) = [U F {r) + U B {r)} N \^) . (5) 

Figure [T] gives a graphical interpretation of this result in terms of a binary tree. It is 
useful to write the expansion of the product on the right side as 

N 

MNt)) = Sn-u,u |^o) (6) 

n=0 

where S m>n represents a sum containing i n+ ^ n ) different terms each comprising n factors 
of U F (j) and m factors of U B (r), and where Q) = k\/[(k — is the binomial 

coefficient. S m ^ n is defined by the recursive relation 

U F (r)U B (kr) (7) 

k=0 
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Fig. 1 Binary tree representation of the generation of the state \^(Nt)) from the origin state |^o) according 
to (J5]l with N = 4. States are represented as nodes (solid discs) and unitary evolution by links (arrows) 
between them. The root node (at the top) represents the origin state |?/>o} and the leaf nodes (on the bottom 
row) represent components of the state |<f r (4r)) where, for example, the node labeled A represents the state 
A.olV'o) = U%\i> ), B represents the state S 3 ,ihM = [U%U F + U 2 B U F U B + U B U F U 2 B + U F U B ]\ip ) and C 
represents So, i[4>o) = U F \ipo). The time axis labeled t at the bottom of the figure represents clock time which 
increases from left to right. Clock time is measured by clock devices whose Hamiltonian H fcl ^ is T invariant. 

with S m fi = UB(mr) and gives, for example, S m) i = J2k=o^B[(fn — A;)t]?7f(t)C/b(A;t). 

The expression (0| J S , jv- n ,n|V'o) represents the evolution of the universe from |^ ) to 
\<f>) over a set of (^) paths through time, where each path comprises a total of n steps 
in the forwards direction and N — n steps in the backwards direction. The set of paths 
includes all possible orderings of the forwards and backwards steps. This set of paths 
is the focus of the remaining analysis. 

First consider the size of the time steps r. To explore the consequences of the limit 
t — > of infinitely-small time steps for a fixed total time t tot , we set r = t tot /N 
and consider (jSJ) for increasing values of N. We can write 2~ n [Uf(t) + Ub(t)] n = 
{exp[— ioiHp — Hb)t\ + 0(t 2 )} n which becomes exp[—i^(Hp — HB)t to t] as N — > oo. In 
this limit, the universe evolves according to the Hamiltonian \{Hp — Hb)- This means 
that isolated subsystems of the universe that obey T invariance, such as our model of a 
clock, would not exhibit any evolution. As clocks are taken to measure time intervals, 
the universe as a whole would not exhibit dynamics in the conventional sense in this 
t — > limit. The lack of dynamics is an unphysical result if our system is to model 
the visible universe. This suggests that the time interval r should be a small non-zero 
number for a universe-like system. Accordingly, in the following we set the value of r to 
be the smallest physically-reasonable time interval, the Planck time, i.e. r pa 5 x 10~ 44 s. 

We now derive a more manageable expression for S m ^ n . The first step is to reorder 
the expression so that all the Ub factors are to the left of the Uf factors using the 
Zassenhaus formula [26]. This yields 

S m , n = U B (mr)U F (nr) £ • ■ ■ £ £ exp [(v + ■■■ + £ + k)r 2 [H F , H B \] exp[(9(r 3 )] (8) 

v=Q e=o k=0 

where there are n summations on the right side and [A, B] is the commutator of A and 



B. |Appendix A| gives the details of the calculation. Next, the analysis is made simpler 
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by expressing the summand in the eigenbasis of the Hermitian operator i[Hp, Hb\- For 
this we use the resolution of the identity given by 1 = / 17(X)dX where the measure 
dX accommodates both continuous and discrete spectra and i7(A) is the projection 
operator that projects onto the state space spanned by the eigenstates of i[H F ,H B ] 
with eigenvalue A. The degeneracy of the eigenvalue A is given by the density function 
p(X) = Tr[i7(A)] where Tr is the trace. To emphasize the degeneracy of A we write the 
identity operator as 

1= f P (X)n(x)dx (9) 



where n(X) = I7(X)/p(X) is an operator with unit trace. Thus, for example, i[Hf,Hb] 
J Xp(X)II(X)dX. Multiplying (JHJ) on the right by ([9]) and ignoring the term of order 



t 3 gives 



where 



S m ,n = U B (mr)U F (nr) J I m>n (A)p(A)i7(A)dA (10) 



4,, b (a) = J2 ■■■ J2 ex p H ( v + • ■ ■ + £ + k ) r2x ■ 

v=0 1=0 k=0 

In the following analysis terms of order r 3 are negligible and are ignored. Performing 



the algebraic manipulations described in Appendix B eventually yields 



T m na{l-expH(n + m-g)r 2 A]} 



There are three important points to be made about the results in Eqs. ( ITU]) and (112]) . 
The first is that if the universe's origin state j^o) is a A = eigenstate of the operator 
i [Hp, H B ] then the integral in (fTOj) can effectively be replaced with I m „(0)p(0)i7(0) 
and the resulting evolution is equivalent to the T invariant case. In order to focus on 
the implications of T violation we assume that the universe's origin state \ipo) has a 
zero overlap with all A = eigenstates, i.e. 

/7(0)hfo>=0. (13) 

We call this the nonzero eigenvalue condition for convenience. We look at the impli- 
cations of relaxing this condition in Section |5j The second point is a useful symmetry 
property of I m ,n{X). Multiplying both the numerator and denominator of the right side 
of (1121) by rifcir [1 — exp(— ikr 2 X)}, canceling like terms and judiciously relabeling in- 
dices yields the symmetry 

WA) = 4, m (A) . (14) 

For the third point, it is straightforward to show for n < m that the zeroes of the 
denominator of I m ,n(X) in (fT2"|) do not result in singularities, because if any factor 
[1 — exp(— iqr 2 X)} in the denominator is zero for A = A , there is a corresponding factor 
{1 — exp[— i{n + m — £)t 2 X}} in the numerator, where £ = n + m — jq and j is a positive 
integer, for which 

1 — exp[— i{n + m — £)r 2 X] 
1 — exp(— iqr 2 X) 

as A — > Xq. The symmetry of I m , n {X) in ( 1T4|) means that the same result holds for the 
complementary case n > m. 
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Fig. 2 Comparison of the func- 
tion |/]v-n,n(A)| (solid curves) 
with the quadratic approxima- 
tion (dashed curves) for various 
values of N and n. All curves 
have been scaled to give a peak 
value of 1. 




3 Interference between paths 

We now analyze the behavior of Sjv-n,n i n (E]) in some detail. Let tp = nr and = 
(N — n)r be the aggregate times that the universe evolves forwards and backwards, 
respectively, where iV is the total number of steps and £ t ot = tp + ts = Nt is the 
corresponding total time. 

Certain features of lN-n,n{ty are more apparent when its modulus |/jv-n,n.(-^)| is 
written as the product of diffraction grating interference functions, i.e. 



'7V-n,n 



(A) | 



j-j-n sin(a g /3 g r 2 A) 



sm((3 q T 2 \) 



where a q = ~(N + 1 — q) and /3 g = |. \lN-n,n(ty\ is a periodic function with a period of 
27r/r 2 and has maxima at A = k2ir/T 2 of n™ =l a q = for integer k. For the times tp = 
nr and ts = (N — n)r to be physically measurable intervals, both n and (N — n) need 
be very much larger than unity. In these cases we can treat |ijv-n,n(A)| as comprising 
narrow peaks at A = k2n/T 2 , for integer k, and being negligible elsewhere. Expanding 
|/jv-n,n(A)| as a power series in A gives the quadratic approximation |/ m>n (A)| ~ (^)[1 — 
j^n(N — n)(A r + l)r 4 A 2 ]. Fig. |2] illustrates the accuracy of the approximation near A = 0. 
The zeroes of the quadratic approximation can be used as a reasonable estimate of the 
width Wn-u n of the peaks. This estimate is given by 



21 



N—n,n 



r 4 n(N - n)(N + 1)] 



To estimate typical values for the eigenvalues A we take neutral kaon evolution as 
a prototypical T-violating process. We first consider the problem for a single kaon for 
which we label the eigenvalues and Hamiltonians with a superscript (1). The eigenvalues 
A^ 1 ) of i[Hp\ H^] in the 2 dimensional subspace spanned by the kaon and anti-kaon 
states can be found using the phenomenological model of Lee and Wolfenstein [T^]. In 
Appendix C| we show that this gives A 1 ^ pa ±10 17 s~ 2 (in units where K — 1) using 



empirical values of Yao et al. 
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Next, we estimate the eigenvalues of the operator i[Hp,H B ] for the universe con- 
taining M kaons. For this we use the fact that the Hamiltonian of a non-interacting 
set of subsystems is simply the sum of the Hamiltonians of the subsystems. Assuming 
that the collection of kaons in the universe do not interact we can approximate each 
eigenvalue A of the operator i [H F , H B ] as a sum A = J2j M where A^ is an eigenvalue 

of i[Hp\ H B ^] for the jth kaon. The values of A range from — MjA 1 - 1 ^ to MIA 1 - 1 ^. The 
degeneracy of the eigenvalue A = A; | A*- 1 ^ | for integer k is given by the density function 
p(A) = M\/(\M - k)\(\M + k)\. Note that 2" M p(A) is approximately a Gaussian dis- 
tribution over A with a mean of zero and standard deviation of \$d — \\fM\\^\ and 
so we shall take the typical values of A to range over —Asd, • ■ • , Asd- Setting M to be a 
fraction, / < 1, of the average total number of particles in the universe, say 10 80 , gives 
M « /10 80 and so A SD ~ v 7 /^ 57 s~ 2 . 

There are two important limiting cases to consider that depend on the relative sizes of 
the width parameters Wjv-n.n and Asd of lN~ n ,n(^) and p(A), respectively, as illustrated 
in Fig. [3j The first is the regime given by Wjv-n.n 3> Asd, i-e. relatively broad peaks 
of ^v-n,n(A). In this limit 7jv_ n)Tl (A) can be treated as approximately constant over the 
range |A| < Asd and so the integral in ([10]) is well approximated by (^)l where 1 is the 
identity operator. This implies that there is negligible destructive interference between 
the paths represented by (^ISjv-n.nl^o)- The second case is the regime W^v_ n ,n ^ ^sd, 
i.e. relatively narrow peaks of lN-n,nW- In the extreme version of this, 7jv- n|Tl (A) is 
zero except in the neighborhood of A = 0. The integral in (flOi) therefore becomes 
(n)II(Q) oc 71(0), where 77(0) is the projection operator which projects onto the zero 
eigenvalue subspace, and so SN- n ,n = Ub[(N — n)r]L r ir(nr)(^ r )77(0). But according to 
the nonzero eigenvalue condition in (1X31) . the origin state of the universe has no overlap 
with any zero-eigenstate state of the commutator and so Sjv-^nlV'o) = 0. In other 
words, the corresponding paths tend to undergo complete destructive interference in 
this regime, and make negligible contribution to the sum on the right side of 

The physical implications of these two limiting cases become apparent on considering 
both forwards and backwards aggregated intervals tp = nr and t B = (N — n)r to be 
very much larger than the Planck time r , i.e. n, (N — n) ^> 1. In this regime the width 
of the peaks of 7jv- n ,n(A) can be approximated as 



\/2A 1 \/24 

W N . n , n = . V [1 + O(-)} 



^Jr^niN - n)N N r l l 2 ^t F t B t tot 

where t t ot = tp + t B = Nt is the total time. Complete destructive interference occurs 
for WN- n ,n < Asd which implies that (tF^tot) 173 > t" _1 ^ 3 (Asd) _2 ^ 3 - For t F ~ t B this 
implies 

*tot ^ ~rjT7\ Ti77 ~ / l73 10 23 s . 
^ i/3 (A S D) i/3 

The destructive interference means that the corresponding set of paths do not contribute 
to the sum on the right side of ([6]); these paths are effectively "pruned" from the set 
possible paths that the universe can take. 

It is the condition tp ~ t B ^> r that is responsible for this destructive interference. So 
to find paths that survive the interference we need to consider times tp t B or t B tp 
for total times t tot = Nt which are at least measurable. For example, consider total 
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Fig. 3 Comparison of the density function p(X) (dotted curve) with the integrand function |ijv- n ,n(A)| (solid 
curve) and its approximation (dashed curve). The curves representing |/jv-n,n(A)| and its approximation have 
been scaled to give a peak value of 1, whereas for clarity, the curve representing p(A) has been scaled to give 
a peak value of 0.5. The standard deviation of the density function is given by Asd = 0.0002/r 2 . The curve 
|ijv-n,«(A)| for TV = 500 and n = 100 has a relatively broad peak whereas for N = 8000 and n = 400 the peak 
is relatively narrow compared to the peak in the density function curve. 



times at the current resolution limit in time metrology [T7j in which case t tot « 10~ 17 s 
and so N ~ 10 27 . Clearly there is no destructive interference for the cases n = or 
n = N (because there is only one path in each case). Also for n ~ 1 and n « N — 1 the 
estimated width of the peaks of lN-n,nW are given by 

\/24 1 

W N „ n , n = ^{l + 0(-)} . (15) 

In these cases Wjv-n,n ~ 10 61 s~ 2 ^> Asd which implies negligible destructive interfer- 
ence. This means that the sum in (EJ) includes a number of paths, however, the paths 
are characterized by one of either tp or ts being negligible compared to the total time 
of 10~ 17 s. 

Finally consider the value of t tot for which there is complete destructive interference 
for the sets of paths represented by (0|Sjv-n , n |^>o) for & h values of n except for n = and 
n = N. In other words there is complete destructive interference for n — 1, 2, . . . (N — 1) 
and so fl§]) is simply \\P(Nr)) = 5^,0 1 V'o) + So.jvIV'o)- We only need to find the value of 
ttot f° r the case n = 1 as this also ensures the destructive interference in the remaining 
cases. Thus we want to know the value of t tot that satisfies Wjv-1,1 <§C Asd- The width 
Wn-i,i is approximated by (fT5l) and so we have y/2A/ (t 2 N) Asd which implies t tot ^> 

y-l/2 10 -13 s In this 

regime, the two paths that survive the destructive interference 
represent exclusively forwards and exclusively backwards evolution. Figure H] shows the 
effect of the destructive interference on the binary tree representation of the evolution. 

In summary, for total times t = Nt > 10~ 17 s (JH]) can be replaced with l&it)) = 
[Euwo&-ii,ii + Z)n«jv'Sjv-n ) n]|'0o) an d the integral in ffTU]) is effectively the identity 
operator. We note that for n w 0, (fTU|) then implies that S'jv-n.n — Ub(Nt) + O(r). 
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Fig. 4 The "pruning" of the bi- 
nary tree in Fig. [T] by destruc- 
tive interference for the case 
Hp ± H B and N > 1. The rela- 
tive weighting of each path seg- 
ment is depicted by the greyness 
of the corresponding arrow with 
black and white representing 
maximal and minimal weight- 
ing, respectively: (a) shows an 
expanded view of the detail near 
the root node and (b) shows the 
whole tree on a much coarser 
scale. 



Similarly, it can be shown using the symmetry property (IT4"|) that S^-n,n = U F (Nr) + 
0{t) for n ~ N . On ignoring the terms of order r we arrive at a key result of this 
paper, the bievolution equation of motion 

\V(t)) = iU B (t) + U F (t)} \i/j ) . (16) 

The term "bievolution" here refers to the dual evolution generated by the two different 
Hamiltonians. The approximations made in deriving this equation become exact in the 
limit t > f-WlO-^s. 

4 Unidirectionality of time and standard quantum theory 

Consider a situation where the Hamiltonians H F and Hb leave distinguishable evidence 
in the state of the universe. For example, imagine that the universe's state at (total) 
time ti is 

1^)) = 10) + 10> 

where \<f>) = U B (ti)\ifj } and |0) = U F (ti)\i/; ) . Let U^fa + a) = V p(a)E/ ^(h) and 
U^(a)\x) = \Xfj) where /i is either F or B, x is either or 9, and a is a suitable 
fixed time interval. Then, according to f)16p 

|!F(ti + a)) = U B (ti + a)|^ ) + U F (h + a)|^ ) 

= Ub^Ub^Hjo) + UpiajUpih)^) 
= U B (a)\(f>) + U F (a)\6) 
= \<Pb) + \0 f ). 

The Hamiltonians leave distinguishable evidence in the state of the universe if they 
generate orthogonal states in the sense that (9b\9 f ) = (<Ab|<Af) = 0. We assume this 
to be the case. Next, let this process occur twice such that U u (a)U fl (a)\x) = \x^ v ) and 
U^(ti + 2a) = U^(a)U^(ti + a) where \i and v are each either F or B. We find, again 
by dlSD, that 
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Fig. 5 The actions of Up(t) 
and f7s(t) on various states of 
the universe in relation to clock 
time t. Diagram (a) represents 
the actions of Uf(t) and Ub(t) 
on an arbitrary state at 
clock time to- Diagrams (b) and 
(c) represent the description of 
the evolution of the universe by 
an observer who finds evidence 
that the Hamiltonian is Hf or 
Hb, respectively. 



(a) 



(b) 



(c) 



\P> 



U B (t) U F {t) 

M 



U F (t 2 ) 



U B (t 2 ) 



H >t 

t +h 



\d') 



+ 2a)) = U B (ti + 2a)|Vo) + U F (h + 2a)|^o) 

= U B {a)U B (t x + a)|^o) + Up^Up^ + a)|^ ) 
= U B {a)\<p B ) + U F {a)\e F ) 

= \<Pb,b) + \0 f ,f) (17) 

which shows that corroborating evidence is left in each term of the superposition. Recall 
that the opposite signs in the exponents of the definitions of Up and U B in Eqs. ([1]) and 
([2]) imply that they represent time evolution in opposing directions of time, as illustrated 
by Fig. [5](a). We conclude that the expressions Up(t)\^) and U B {t)\^) represent the 
universe evolving in opposite directions of time and in each case the evolution leaves 
corroborating evidence of the associated Hamiltonian in the state of the universe. 

To interpret what these results mean, consider an observer within the universe who 
performs measurements to determine the nature of the Hamiltonian. According to ( TT7|) 
the observer would be in a superposition of two states, one representing the observer 
consistently obtaining physical evidence that the Hamiltonian is Hp, and the other 
representing the observer consistently obtaining physical evidence that the Hamiltonian 
is H B . This means that each term on the right side of the bievolution equation (TIB]) 
represents the observer having access to only one version of the Hamiltonian. It follows 
that the observer would describe the universe as evolving according to either 

\$(t)) = U F (t) |^o) (18) 

or 

\V(t)) = U B (t) |Vo) (19) 

depending on whether the observer has evidence that the Hamiltonian is Hp or H B , 
respectively. Each of (ITS]) and (1T9"|) is the solution of the conventional Schrodinger equa- 
tion for the corresponding version of the Hamiltonian. The observer's description would 
therefore be in accord with our own experience of the dynamics of the universe. 

It is worth remarking that there is no inherent preferred direction of time in conven- 
tional quantum mechanics. Indeed, in conventional quantum theory the Hamiltonian 
is the generator of time translations in either direction of time. Here, however, each of 
(TTS]) and ffTS]) is associated with a particular direction of time evolution in the sense that 
Up{r) evolves the state of the universe along the positive— t ("forward") direction and 
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Ub{t) evolves the state of the universe along the negative— t ("backward") direction, as 
illustrated in Figs. EJ^b) and|5^c). This association stems from Hp being the generator 
of time translations exclusively in the positive— t direction and Hb being the generator 
of time translations exclusively in the negative— t direction in accord with Wigner's 
time reversal operation [TJ] as discussed in Section 2. The association has important 
consequences for the order that events occur in time. For example, consider the case 
where the observer finds evidence that the version of the Hamiltonian is Hp as shown in 
Fig. E(b). Time evolution in this case is in the positive— t direction. Imagine that a par- 
ticular atom is undergoing spontaneous radiative decay at t — to and let |e) represent 
the universe with the atom in an excited energy eigenstate. The state \d) = Up{t,2)\e) for 
an appropriate value of £2 would represent the universe with the atom having decayed 
to a lower energy eigenstate along with its emitted radiation, say. The state \p), where 
\e) = Up(ti)\p), could represent a stage in the preparation of the excited atom. The 
order that the preparation stage, excited atom and decay atom appear in time is dic- 
tated by the direction of the time evolution. The converse case where the observer finds 
evidence that the Hamiltonian is Hb is illustrated in Fig.|^c) for which \d'), \e') and \p') 
represent the decayed atom, the excited atom, and the preparation stage, respectively. 
The reversed ordering of the events in this case is due to the opposite direction of time 
evolution that is associated with the Hb version of the Hamiltonian. In both cases the 
observer would find physical evidence of a fixed direction of time evolution and that 
direction would correspond to the version of the Hamiltonian observed. 

We have allowed the universe to take any path through time according to (JSJ). But 
provided there are sufficient T violation processes, destructive interference between 
paths will result in the universe evolving according to the bievolution equation (IT6|) . 
We have shown that in this case an observer within the universe would find physical 
evidence of only one version of the Hamiltonian. As a consequence, the observer would 
describe the universe as evolving according to the Schrodinger equation in the direction 
of time that corresponds to the observed version of the Hamiltonian. We have therefore 
demonstrated the consistency between the bievolution equation and conventional quan- 
tum theory for an observer within the universe and how the unidirectional nature of 
time arises from the perspective of the observer. 

5 Eigenstates of i[H F , H B ] with zero eigenvalues 

We now temporarily relax the nonzero eigenvalue condition in f JT3|) and set the origin 
state to be an eigenstate of the commutator i[Hp, Hb] with eigenvalue zero, i.e. we set 

7T(o)|^o> = |^o> • 

As previously mentioned, the integral in (fTO]) in this case can effectively be replaced 
with J m>n (0)p(0)i7(0) = ( n + m )7T(0) and so ® becomes 

N 

\¥{Nt)) = £ U B [(N - u)t]U f {ut) (?) |^o) 

n=0 

which represents all 2^ possible paths through time. We now reapply the analysis of the 
previous section for the processes Up (a) and U B (a) with t 1 = Nr and a = kr. There 
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is no destructive interference in this case, and so the state of the universe at the total 
time ti + 2a is, according to 



\&(t 1 + 2a)) = [U b {t) + U f (t)] 



2k 



■ N 

E 

.n=0 



U B [(N - n)r]U F (nr) 



N 



which represents 2 N+2k paths through time. These paths include the operator combi- 
nations U F (a)U F (a), U f ((i)Ub{q), Ub{o)U f {o), Ub(o,)Ub((i) (among many others) and 
so any distinguishable evidence of the Hamiltonians left in the state of the universe is 
not necessarily corroborated over repetitions of the process. This is not in accord with 
present experimental evidence. Moreover, consider the origin state of the most general 
form 



where /7(0)|^) = and /7(0)|V>H) 
is replaced with 



\ipo). For this case the bievolution equation (JT6 



Mi)) = [Unit) + U F (t)} + E U B [(N - n)r]U F {nr) (f ) \4) 



(20) 



n=0 



where t = Nr. However, as we have already seen, experimental evidence suggests that 
the universe is not following the corresponding set of paths represented by the last term 
on the right side. In order for our system to provide a physical model of the observable 
universe, it is therefore sufficient to only consider origin states that satisfy the nonzero 
eigenvalue condition in ffT3l) . We assume this to be the case in the remainder of this 
work. 



6 Schrodinger's equation for bievolution 

We can construct the differential form of the bievolution equation of motion as fol- 
lows. Increasing t in (fT6|) by a relatively small time interval St yields \\P(t + St)) = 
[6 r ^(5t)L r s(t) + U F (St)U F (t)} I^q) an d so the rate of change of the state is given by 

= [iH B U B (t) - iH F U F {t)\ |^o) + 0(5t) 

where 8\\P(t)) = \\P(t + St)) — \\P(t)). Taking the limit St — > r and ignoring a term of 
order r gives another key result, the Schrddinger equation for bievolution 

d\#(t)) _ d\Mt)) d\Mt)} 

dt dt dt 

where (d/dt)\^(t)) = -iH^(t)) and \^(t)) = U„(t)\^ ) for fx = F or B. 
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7 Discussion 

To gauge the full impact of T violation processes, it is instructive to compare the above 
analysis with a universe which obeys T invariance. In this case Hp = Hb = H and it 
can be shown that 



which is never zero and so all possible paths are included in fl6]). The direction of time 
is ambiguous as there is no physical evidence of any kind to single out one direction 
over the other. The fact that this ambiguity is removed for a universe with T violation 
processes leads directly to the proposition that T violation processes are responsible for 
the phenomenological unidirectionality of time that we observe in the universe. 

Just as important is the implicit assertion here that in the absence of a unidirectional 
nature, time is randomly directed at any given instant. We allowed for the possibility 
of a random direction when we added amplitudes in fl3]). In hindsight, the analysis of 
the preceding sections can be seen to be a study of the dichotomy between randomly 
directed and unidirectional time evolution. In this context the important question about 
the unidirectional nature of time is not how one direction of time is chosen over the 
other, but rather how each direction is maintained consistently from one instant to the 
next. This is one of the underlying questions that is being addressed in this article, with 
T violation providing a possible mechanism as the answer. 

We can also elaborate on what the time reversal operation means in the presence 
of T violating processes. The bievolution equation ffT6]) shows that as the value of t 
increases, the states Up(t)\ipo) and |^>o) m each branch of the superposition trace 
out respective trajectories in the state space. Decreasing values of t then simply "back 
tracks" the states along their respective trajectories in a consistent manner, essentially 
tracking back over what might be called their "histories". In contrast, the time reversal 
operation essentially interchanges Up(t) and Us{t) in ffTEj) . This implies that the time 
reversal operation refers to switching between the branches of the bievolution equation 
as opposed to tracking back over histories. 

These results significantly elevate the importance of the matter-antimatter arrow in 
the study of the nature of time. The underlying T violation is based on the empiri- 
cal evidence that different versions of the Hamiltonian operate for different directions 
of time evolution. Moreover, as discussed in Section HJ each direction of time can be 
uniquely identified by physical evidence left by the Hamiltonian. This is quite different 
to other arrows of time. In particular, the thermodynamic arrow does not distinguish 
between the different directions of time evolution in this way because entropy increases 
in the direction of time evolution regardless of the direction. This suggests that T viola- 
tion plays a far more significant role in determining the nature of time than previously 
imagined. 

We should mention that the physical implications of T violation described here 
depend the magnitudes of the eigenvalues A. Our estimate of these eigenvalues depends 
on the number of particles involved in T violations processes and this has been fixed 
only at /TO 80 where / < 1 is a parameter of undetermined value. However, our analysis 
does not depend critically on the value of / and, indeed, a large range of its values 
support the key results. Also, our analysis has focused on relatively large time scales 
in order to draw out the main consequences of T violation for the unidirectionality of 
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time. In doing so we have ignored terms of order r, the Planck time. These terms imply 
that the directionality of time has finer details for relatively small time intervals. The 
elaboration of these details is, however, beyond the scope of this article. 

There are potentially important experimental implications of the analysis presented 
here. A sufficiently small value of / would produce incomplete destructive interference 
and a deviation from the dynamics described by the bievolution equation ( 116]) . This 
would be manifest in a lower bound to the accuracy of time metrology. A more extreme 
case occurs in the very early universe during the inflation and radiation-dominated 
periods [TS1[TH]. Baryogenesis would not have yet occurred and CP and T violation 
would have been relatively rare events. This implies that the value of / would have been 
negligible during these periods and so the direction of time evolution would have been 
uncertain. This has potentially important ramifications for cosmological models and 
their testing against observational data of the cosmic microwave background radiation 

The analysis presented here may also have important consequences for quantum 
gravity. Consider again the situation where Hp = Hb = H . Instead of expanding ([5]) 
into terms of the kind in (12 lj) we write Up + Ub = 2 cos(Ht) and so 

\V(Nt)) = 2 N cos N (Ht) |^o) • (22) 

Let the eigenbasis of H be given by 

H\E n ,\) = E n \E n ,X) , (23) 

where E n are the energy eigenvalues with Eq = and A indexes different eigenstates in 
degenerate manifolds. Expanding \ipo) in this basis gives 

\#{Nt)) = 2 N Y,c n ,xcos N (E n r)\E n , A) (24) 

where c nt \ = {E n ,\\ip ). Let |^o) be a state of bounded energy to the extent that 
c n ,A = for \E n \ > 1/r. Here 1/r is the Planck energy in units where h = 1. This 
bound would be reasonable for a very young and very small universe. In this case, 
provided c 0; a 7^ for some values of A, the state \\P(Nt)) becomes proportional to 
J2x c o,x\Eo, A) as N — > oo. Hence in the limit as N — > oo 

H\V(Nt)} = . (25) 

According to (J5J), \&(Nt)) is a superposition of states each of which represents a net time 
ranging from —Nr to +Nt. It therefore represents the whole history of the universe. 
( 125]) is analogous to the topological invariance condition discussed by Misner [20] in 
relation to the quantization of gravity. Topological invariance is needed in quantum 
gravity to maintain the invariance of physical quantities to transformations of spacetime 
manifolds. To make a direct connection with Misner's work, we let H represent the 
Hamiltonian H(x) for gravitational and matter fields defined on a four dimensional 
spacetime manifold x, {<r fc } represent a set of spacelike hypersurfaces indexed by integer 
k and separated by timelike intervals of r, Up{tit)[iPq) represent the state functional 
if; an and U B{nr)\ipQ) represent the state functional ^ -_ 7l for positive integer n, where 
ip a is a functional of the gravitational and matter fields at the hypersurface a. In this 
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representation \\P(Nt)) represents a state functional which is a superposition of ijj an 
for n = —N to N. (1231) then represents the result that topological invariant states are 
attractors of time symmetric evolution. The zero energy eigenstate implies that the 
positive energy of matter fields is balanced by the negative energy of the gravitational 
potential energy [2T1I221I23] . 

This means that (I25p corresponds to the Hamiltonian constraint of the Wheeler- 
DeWitt equation [SUHS]. The Wheeler-DeWitt equation represents the entire history 
of the universe with no bias towards either direction of time. Here time symmetric evo- 
lution is unbiased in the same way in the sense that it represents evolution without 
a fixed direction of time and the state in (1251) represent the entire history of the uni- 
verse. These common features suggest that there is an elemental relationship between 
the Wheeler-DeWitt description and the present analysis. The present analysis may 
therefore offer a means of incorporating the two versions of the matter Hamiltonian in 
a Wheeler-DeWitt description of the universe for situations where Hp ^ Hb- Moreover, 
as the present analysis shows how the phenomenological unidirectionality of time arises 
when Hp ^ Hb-, it may also shed fresh light on the problem of how time is established 
in the timeless Wheeler-DeWitt description. 

In conclusion, the physical significance of the relatively weak and rare T violation 
processes exhibited by mesons has been a puzzle for more than four decades. Despite 
representing a fundamental time asymmetry, these processes have been regarded as 
having an insignificant affect on the physical nature of time. This has been due in part to 
the lack of a formalism for accommodating their time asymmetry in a single dynamical 
law. We have removed this obstacle by deriving, from a first principles analysis based 
on Feynman's sum over histories, a general method for incorporating both versions of 
the Hamiltonian for a T violating process, one for forwards and the other for backwards 
evolution, in a single dynamical equation of motion. Moreover we have shown that 
these processes can affect the time evolution of the universe on a grand scale. Indeed, 
we have shown that T violation provides, at least in principle, a physical mechanism 
for the phenomenological unidirectionality of time. Finally, the analysis presented here 
may also have important ramifications for time metrology, quantum gravity and the 
early development of the universe. At the very least it gives a clue to the quantum 
nature of time itself. 
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Appendix A: Reordering factors of S m ,n 
The operator S m , n is defined recursively in (J7|) as 



Til 





with S, 



m,0 — 



Usirnr) and so 



m m 



m,l — 



S m _ kfi Up(T)U B (kT) = £ U B [(m - k)T]U F (T)U B (kT) . 



(A.2) 
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The Zassenhaus formula [26] 

e Ae e Be e -±e 2 [A,B] e O(e3) = & Be & Ae e ~^[B,A] ^(e*) 

(which is related to the Baker-Campbell-Hausdorff formula) where A and B are arbi- 
trary operators and e is a small parameter, can be written as 

e Ae e Be = e Be e Ae e e*[A,B] e O(e>) 

which shows that 

e -ijrH F e ikrH B _ ^krH B e -ijrH F e jkr 2 [H F ,H B ] e O(r 3 ) g\ 

In terms of the operators U F {j) = e~ lTHF and Ub{t) = e trHB this result reads 

U F (jr)U B (kr) = U B (kT)U F (jT)e jkT2[HF ' HB] e° {TS) (A.4) 

or equivalently 

U B (kr)U F (jr) = U F (jT)U B (kT)e- jkr2[HF ' HB] e°( Ti) . (A.5) 

Hence 



Sm,n ^ ] S m — k,n— 1 

U F (r)U B (kr) 



k=0 



,kr 2 lH F ,H B }n(T 3 ) 



= J2 S m-k,n~iU B (kT)U F (r)e k7 

k=0 
m m—k 

= E E ^- fc - J ,n-2f/F(r)t/B(jr)^(fcr)t/ F (r)e^ 2 ^^e ^) 

k=0 j=0 
m m—k 

= E E S m . k . hn . 2 U F {r)U B [{j + k)r}U F (r)e k ^ HF ' H ^e°^ 

k=0 j=0 
m m—k 

= E E S m _ k _ jin _ 2 U F (2r)U B [(j + k)T}e^ 2 l*^ B ] e o(T*) 

k=0 j=0 
m m—k 

= E E S'm-fc-i.n^BlO' + k)r]U F {2r)e^ k ^ H ^e°^ . (A.6) 

fc=0 j=0 

Setting i = j + k and rearranging the order the terms are summed to eliminate j then 
yields 

m i 

Sm, n = E E 5 m _ £ ,„_ 2 f/ B (£r)f/ F (2r)e^ 2 ^'^]e ^) 

m 



«=0 k=0 

771 771 — £ 

E E W3^Wf/B[(i + r)r]U F {2r)e^ H ^ E e^'^e ^ 

f=0 r=0 fc=0 

771 771 — 

E E ^-r,n-3%[(^ + r)T]U F (3T)e^ 2 ^lH F ,H B] £ ^[H^]^)^ 



fcO r=0 k=0 
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Rearranging as before but with j = £ + r and eliminating r yields 

m j I 

S m ,n = E E ^(jr)I/A)« (WlM E e^^e^ 3 ) 

j=0£=0 fc=0 
m j £ 

= E ^-^-3f/B(ir)^(3r)e^^'^] E E e^^leW . (A.8) 
j=o e=ok=o 



Continuing in this way gives, after w uses of the recursive relation (1 A. 1 [) 



S m , n = E S m ^ n . w U B {vr)U F {wr)e^ H - H ^ E " ' E E e^ +e+k ^ H ^e°(^.9) 

v=0 r=0 £=0 k=0 

in which there are w summations. Letting w = n then yields 

7TL V S £ 

S m , n = E S m ^ U B (vT)U F (nT)e^,H B] e (r+-+« + *)T-[^] e o^) 

j;=0 r=0 £=0 k=0 



E f/ B (mr)f/ F (nr)^ 2 ^' fffl l E " ' "E E e^-^+^^'^e ^) 

,(v+r+-+£+k)r 2 [H F ,H B ]n(T 3 ) 



v=0 r=0 £=0 fc=0 

Tfi V s £ 

U B (mr)U F (nr) E E ' " E E e^ +r+ - +e+k ^ H ^e°^ . (A.10) 

u=0 r=0 1=0 k=0 



Hence we arrive at (El) 



S m ,n = U B (mr)U F (nr) E " ' E E e^^^^e ^ (A.ll) 

u=0 1=0 k=0 

in which there are n summations. 
Appendix B: Simplifying J m , n (A) 

Consider the following manipulations of nested summations of the kind 

w j £ s n 

E E E ■ " E E e i{ - j+f - +m+ '" +n+k)e . (b.i) 

j=0 1=0 m=0 n=0 k=0 

We call the summations "nested" because the upper limit n of the last summation over 
k is the index of the second-last summation and so on. The nested property means that 
summations cannot be evaluated independently of each other. However it is possible to 
reduce the number of nested summations by performing certain operations which we 
now describe. 

First, by changing the order in which the indices j and i are summed, we find 

w j l s n w w i s n 

'''EE e i(j+£+m+ "' +n+fc ' )£ ' = e E E '''EE e i ^ +i+m+ '" +n+k ^ (b 2) 

j=0 1=0 m=0 n=0 k=0 1=0 j=£ m=0 n=0 k=0 
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next, by cyclically interchanging the indices in the order j — > k — > n — > s —>■••• — > 
m — >■ I — > j on the right-hand side, we get 

w j £ s n w w j £ s 

y y y y y e i(i+^+m+---+n+fc)e = E E E E ' ' ' E e »(*H-i+<+m+-+n)0 3^ 

J= 0f=0m=0 n=0fc=0 j=0 )t=j 1=0 m=0 n=0 

and finally, bringing the sum over k to the extreme right on the right-hand side gives 

w j I s n w j i s w 

E E E • • • E E e i(j+t+m+-+n+k)e _ y y y . . . y y e iQ-+M-m+-+n+fc)0 ^ 

j=0<=0m=0 n=0fc=0 j=0f=0m=0 n=0 k=j 

We can abbreviate this general summation property as 

w t s n w t s w 

y . . . y y y e i(j+-+s+n+fc)e _ y . . . y y y e i(i+-+s+n+fc)e _ 

j=0 s=0n=0fc=0 j=0 s=0n=0k=j 

Consider the product 

S TL S TL S TL 

(jo + 1) e e ei{n+k)d = e i9 EE el{n+k)e + EE ei{n+k)e ■ ( B - 6 ) 

n=0 k=0 n=0 k=0 n=0 fc=0 

Using flB.5j) for the two summations in the first term on the right-hand side gives 

s n s s s n 

(e ie + l) E E e i(n+k)e = e* 9 E E e i{n+k)e + EE e i{n+k)e 

n=0 k=0 n=0 k=n n=0 k=0 

s s s n 

= y V e i ( n + k + 1 ) e + E E e i(n+k)e 

n=0 k=n n=0 k=0 

s s+1 s n s s+1 

= E E e i(n+k)e + y y e i(n+fc)e _ y y e i(n+fc)0 ^ ^ 

n=0 k=n+l n=0 k=0 n=0 k=0 

The two nested summations on the left-hand side have been reduced to two un-nested 
summations on the right-hand side. Similarly, using (1B.5|) again as well as (1B.7j) we find 

t s n t s n 

D i{s+n+k)8 



( e fifl +e * +1 ) J2J2H e^ +n+ W = e i2e E E E 

" =0 

■( e,9 + 1 )EEE e!(s+n+fc)9 



e 

s =0 n=0 A;=0 

t s n 



s=0 n=0fc=0 
t s t t s s+1 

_ y y y e i(s+n+k+2)e + E E E e i( - ,s+n+fc ^ 

s=0 n=0 fc=s s=0n=0t=0 
t s t+2 t s s+1 

_ y y y & i{s+n+k)e + y y y e i(s+n+k)9 

s=0n=0k=s+2 s=0n.=0fc=0 
t s t+2 

= E E E e i{s+n+k)e . (B.8) 

s=0 n=0 fc=0 

Here the t/iree nested summations on the left-hand side have been reduced to two nested 
summations and one un-nested summation on the right-hand side. This implies that, 
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for example, p + 1 nested summations can be reduced to p nested and one un-nested 
summations as follows: 



p w t s n w t s w+p 

' Y eir6 ) Y ■ " Y Y Y e i(i+ '" +s+ri+fc)e = ( y ■■ ■ y Y e i(i+ " +s+n)e ) y < 

>=0 j=0 s=0n=0fc=0 j=0 s=0n=0 fc=0 



p+1 nested sums p nested sums (B-9) 

Moreover, by repeating this process once more, the p nested summations of the right- 
hand side can be further reduced to p — 1 nested and one un-nested summations. 
Continuing in this way eventually leads to the following result 

p q w t s n p—1 q w t s w+p 

II ( Y eire ) E"'EEE e in+-+s+n+k)e = y\[y eire ) ( Y" ~Y Y e i(j+ "' +s+n)e ) y eike 

5=0 r=0 j=0 s=0n=0t=0 <?=0 r=0 j=0 s=0n=0 k=0 

p+1 nested sums p nested sums 

w w+p— 1 w+p 

= (EO"-( Y e™°)(Ye lke 

j=0 n=0 fc=0 



Thus 



p+1 un-nested sums 

p w+q 

n(E^)- (b.io) 

q=0 j=0 



w t s n Y[ p rA T w+q e^ e 

e i(j+-+s+ n +k)e = iig=oy^=o j , B , 

^=0 s=0n=0k=0 1 L=0 I <Lr=0 e 



p+1 nested sums 

Evaluating the two geometric series on the right-hand side then gives 



w t s n 

V ryy e i(j+-+s+n+fc)e = 11 g=° v )_ = n i =1 v 



p+1 nested sums (B.12) 

Adjusting the number of nested summations of the left-hand side then gives the useful 
result that 

w t s n n P _i ( 1 — e^™" 1 " 9 ) 6 ^ 

^ . . . Y Y e iU+ '" +s+n+h)e = - - . (B.13) 



j=0 s=0 n=0 fc=0 



p nested sums 



We can now use this result to simplify the n nested summations in flTT]) as follows. First 
we note that 

m si yt n _ (l — e ~i{m+q)T 2 \\ 

WA) = Y---YY e-^-+^ = i^ii — ; (B14) 
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and then by redefining the index q in the numerator we arrive at 

J~jra— 1 ^ g— i(n+m— <?)t 2 A n 



i(A) 



n" i (1 - e-n-^ 



(B.15) 



which is ffT2j) of the paper. 

Appendix C: Estimating eigenvalues 

We estimate the eigenvalues of the commutator [Hp \ on the state space of a kaon 
field mode that represents a single particle. However, to avoid unnecessary clutter, we 
will omit the superscript "(1)" on the single particle Hamiltonians in the remainder of 
this section. 

It is convenient to invoke the CPT theorem and use the CP operator instead of the 
T operator as follows 



H B = THpT- 1 = (CP)- 1 H F (CP) . 



(C.l) 



Consider the matrix representation in the kaon state space where the vectors (1,0) and 
(0, 1) represent the kaon and anti-kaon states \K°) and \K°), respectively. In this basis 
the general form of the CP operator is 



CP 



e 



ill 







(C.2) 



where t] and £ are possible phase angles. By redefining the complex phase of the kaon 
state where \K°) is multiplied by e^ v ~^^ 2 , the CP operator becomes 



CP 



e M 
e id 



jo 



1 

1 



(C.3) 



where 9 = -(77 + £). (This is a more useful transformation for us here compared to the 
usual one, which multiplies \K°) by e - ^.) Let be the matrix elements of Hp in 
the kaon state basis where Mn = (K°\Hp\K ) etc. To find the corresponding matrix 
elements of Hp we use fIC.ll) . that is 

H B = (CPy l H F (CP) 





"o I 




"M n 


M 12 ~ 




"0 l" 




' M 22 


M 21 ~ 




1 




M 21 


M 22 




1 




M 12 


M n 



Hence we immediately find 
HpHs = 
H B H F = 



M n M 22 + (M 12 ) 2 M n (M 21 + M 12 ) 
M 22 (M 21 + M 12 ) M U M 22 + (M 21 ) 2 

M n M 22 + (M 21 ) 2 M 22 (M 21 + M l2 ) 
M U (M 21 + M l2 ) M n M 22 + (M l2 f 
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and so the commutator is 
[Hp, Hb] = HpHs — HbHf 

which is in the form 



(M 12 ) 2 - (M 21 ) 2 (M u - M 22 )(M 21 + M 12 
■(Mn - M 22 )(M 21 + M 12 ) -[(Mi 2 ) 2 - (M 21 ) 2 ] 

(C.4) 



[H Fl H 



A B 
-B -A 



(C.5) 



The eigenvalues of the last matrix are ±V A 2 — B 2 . Note that A is imaginary and B is 
real, and so the eigenvalues are imaginary. We now determine empirical values for the 
variables M isj in flC.4[) . 

The phenomenological model given in Lee and Wolfenstein [TH] describes the evolu- 
tion of the amplitudes a(i) and b(t) for the single-particle kaon field to be in the states 
| A' ) and \K°), respectively, by the differential equation 



dij){t) 



[r + iM)ii>{t) 



dt h 

where ip = The matrix on the right-hand side is defined by 

r + iM = D + i{E x a x + E 2 a 2 + A 3 a 3 ) 
where the matrices r and M are Hermitian and 





"o r 




"0 -i 




"1 " 


0"! = 


1 


0"2 = 


% 


0-3 = 


-1 



(C.6) 



(C.7) 



(C.J 



The analysis below shows that the matrix M in (1C6|) represents a CP violating Hamil- 
tonian. For our purposes, it provides a sufficient estimate of the T violating Hamiltonian 
for a neutral kaon. Accordingly, we use the empirical values of the elements of M as 
estimates of the corresponding matrix elements Mij of the Hamiltonian Hp that ap- 
pears in (1C.4j) . We should also mention that we could use M as an estimate of the 
matrix elements of Hb instead, however the only difference this would make is that the 
matrix elements of the commutator [Hp,Hs] would be multiplied by —1, and this has 
no physical consequences because the two eigenvalues are the negatives of each other. 

If the CPT theorem holds, which we assume to be the case, then E\ = Acos0, 
E 2 = EsiiKp and A3 = 0. Here is a complex parameter [16] which is related to the 
CP violation parameter e by e = (1 — e 1 ^)/ (1 + e 1 ^), or equivalently, 



j4 



1 + e 



(C.9) 



The values of D and A comprise the decay rates 7$ and masses m, of the short and 
long-lived components, |A°) and |A°) respectively, of |A°) as follows 



D 

1E 



I (71 + + :H m i + m 2) 

|(Ti - 72) + \i{m x - m 2 ) . 



(CIO) 
(C.ll) 
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Here 7 $ and rrii are the energies 7 j = and rrii = m^c 2 , respectively, where 7^ and rri! i 
are the actual decay rates and masses, respectively. Also note [27] that |e| ~ 2.3 x 10~ 3 
and arg(e) ~ 45° and so from (1C.9j) we find e 1 ^ = 1 — 2e + 0(e 2 ). We can therefore make 
the following approximations 

e^ = l-v / 2|e|(l + + O(e 2 ) , 
e ir = l /^4>y = 1 + 72| e |(l - i) + 0(e 2 ) , 
e -M = i/( e ^) = 1 + V2je](l + z) + 0{e 2 ) , 
e -ir = fjty = ! _ y/2\ e \(l - i) + 0(e 2 ) . 

This allows us to write down the following useful results 

= 1 - zV2|e| + 0(e 2 ) , 
= 1 + iV2\e\ + 0(e 2 ) , 



i(e i *-e^) = -v / 2|e| + C(e 2 ) 



-10* 



V2\e\ + 0(e' 



From Eqs. (ICTfjl . (El, (jCTty and (lUTTj) we find 
r + iM = D + z'(£ cos (Ti + £ sin cr 2 ) 

|(7i + 72) + + m 2 ) |( 7 i - 72) + \i{rn\ - m 2 ) 



-itj> 



|(7i -72) + -m 2 ) 

The Hermitian transpose of this equation is 
r-iM = 



7(71 + 72) + \i{m x + m 2 ) 



(C.12) 



- 72) - \%{m x - m 2 ) 
(71 - 72) - \i{m\ - m 2 ) e l<t> * \^f x + 72) - \i{m x + m 2 ) 



l(7i +72) - \i{m x +m 2 ) 



(C.13) 



By subtracting (IC.13P from (10. 1 2j) and dividing the result by 2i we get 



M 



|(mi + m 2 ) 
- 7 2 )V2|e| + i(mi - m 2 )(l - z>/5|e|) 

f C(e 2 ) . 



- 72)V2|e| + |(m! - m 2 )(l + z>/2|e|)' 



|(mi + m 2 ) 



(C.14) 



This shows that M xx = M 22 and so S = in Eqs. (1C.4j) and (1C.5I) . The eigenvalues of 
the commutator [Hp, Hb] are therefore simply ±A where 



A = M 2 2 - M 2 2l 



zi( 7l - 72 )V2|e| + \{m x - m 2 )(l + iV2|e|)] 2 

z±( 7l - 7 2 )V2|e| + |(mi - m 2 )(l - zV2|e| 
-iz(7i - 72)V2|e||(mi - m 2 ) [2(1 + iy/2\e\) + 2(1 - iV2|e|) 



0{e 2 



+i(mi - m 2 ) 2 [(1 + zV2|e|) 2 - (1 - zv^lel) 2 ] + 0(e s 
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= i[ - (71 - T2)V2|e||(mi - m 2 ) + ( mi - m 2 )V2|e|] + 0(e 2 ) 
= zV2|e| [( 7l - 7 2 )K m 2 - mi) + (mi - m 2 ) 2 ] + C(e 2 ) 

= zV2|e| [±Z\ 7 Z\m + {Am) 2 ] + C(e 2 ) (C.15) 

where Z\m = (m 2 — mi) and ZV7 = (71 — 72). 
Finally, the empirical values [27] 

Am « 0.56 x lO 10 ^" 1 , 

Ay « 1.1 x lO 10 ^" 1 , 

show that the eigenvalue A in (1C.15j) can be approximated as 

A m H0 17 h 2 s~ 2 . (C.16) 

This is the value used in the paper for \^ where A^ = ±iA. Note that in the main text, 
we use units in which h — 1. This is equivalent to replacing the Hamiltonians Hp and Hb 
in this section with Hp/h and Hb/H, respectively. In that case the commutator in (10. 5|) 
has the form [Hp, Hb]/Ti 2 which has a dimension of (time) 2 . Hence A^ ~ ±10 17 s" 2 . 
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